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We analyze spontaneous parametric down-conversion in various experimentally feasible ID 
quadratic nonlinear waveguide arrays, with emphasis on the relationship between the lattice’s topo¬ 
logical invariants and the biphoton correlations. Nontrivial topology results in a nontrivial “winding” 
of the array’s Bloch waves, which introduces additional selection rules for the generation of bipho¬ 
tons. These selection rules are in addition to, and independent of existing control using the pump 
beam’s spatial profile and phase matching conditions. In finite lattices, nontrivial topology produces 
single photon edge modes, resulting in “hybrid” biphoton edge modes, with one photon localized at 
the edge and the other propagating into the bulk. When the single photon band gap is sufficiently 
large, these hybrid biphoton modes reside in a band gap of the bulk biphoton Bloch wave spectrum. 
Numerical simulations support our analytical results. 
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I. INTRODUCTION 

Spontaneous parametric down-conversion (SPDC) is 
an important way to generate pairs of photons exhibit¬ 
ing quantum correlations, with applications ranging from 
fundamental tests of quantum theory (Bell tests) to 
quantum cryptography and information processing |l|-l3j]. 
Genuinely quantum behaviour and scalability require a 
high fidelity of the photon pairs, which is limited if they 
are generated and shaped by bulk optical components. 
Hence there is currently strong interest in implementing 
SPDC in integrated optical devices M- 

SPDC in nonlinear waveguide arrays has been pro¬ 
posed as a tool to tailor biphoton quantum correlations 
in integrated optics (Ml- It offers many advantages 
over bulk components; biphotons are generated directly 
in the device, so there are no input coupling losses. The 
biphoton spectrum and correlations can also be read¬ 
ily controlled via the pump beam’s spatial profile and 
the diffraction or quantum walk of the generated pho¬ 
ton pairs through the array 0®. This concept was 
recently demonstrated in experiments in lithium niobate 
waveguide arrays ®®. 

So far however, SPDC was only studied in homo¬ 
geneous waveguide arrays, or arrays with a single de¬ 
fect H71 |. In both cases only a single Bloch band 

is relevant. It is interesting therefore to explore the op¬ 
portunities for controlling SPDC and biphoton correla¬ 
tions offered by modulated arrays with multiple Bloch 
bands (2pj |. 

Given the endless possibilities in designing modulated 
waveguide arrays, it is useful to group them into differ¬ 
ent classes sharing similar properties. One way to do 
this is using topological invariants , which has lead to 
a wide range of breakthroughs both in the fundamen¬ 
tal band theory of solids, to devices using “topologically 
protected” surface states that are robust against disor¬ 
der (2ll-[23l| . These ideas are now attracting interest in op¬ 


tics, and several photonic analogues of these topological 
condensed matter systems were recently demonstrated in 
experiments ®. 

Here we explore how lattice topology can provide an 
additional tool to control correlations of biphotons gen¬ 
erated in quadratic nonlinear waveguide arrays. The 
basic idea is that nontrivial topology is associated with 
a nontrivial “winding” of the lattice’s Bloch waves, and 
this winding can lead to selection rules controlling which 
biphoton modes are strongly excited. In principle, the 
mode winding can be completely independent from the 
array’s dispersion relation (phase matching conditions), 
so it offers another degree of freedom to control biphoton 
correlations in integrated optics. This topology is inher¬ 
ently robust against fabrication disorder. In some cases 
lattices with nontrivial topology also host protected edge 
modes, which produce bands of “hybrid” biphoton states 
exhibiting entanglement between localized and propagat¬ 
ing modes. We demonstrate the feasibility of these ideas 
by carrying out numerical simulations of various one di¬ 
mensional (ID) topological lattices under the tight bind¬ 
ing approximation. The model and parameter regimes 
are accessible in current state of the art experiments. 

In Sec. El we review the theory of SPDC in waveguide 
arrays, generalising to multi-band (modulated) arrays. 
In Sec. m we make some general statements on the role 
of topology in two band, ID models. Following this, as a 
concrete example we consider in detail the Su-Schrieffer- 
Heeger (SSH) model in Sec. HVl which illustrates main 
features and is a practical, experimentally realisable ex¬ 
ample. Sec. E| compares these results against a binary 
lattice where the waveguide depths are modulated, which 
is an example of a “nontopological” model because it lacks 
the required symmetry. We conclude in Sec. I VI I with a 
summary and discussion of future directions. 
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FIG. 1. (Color online) Schematic of quadratic nonlinear 
waveguide array. Pump beam generates photon pairs that 
diffract to neighboring waveguides. Quantum correlations of 
the photon pairs can be measured via coincidences of two 
single photon detectors. 


II. SETUP AND OBSERVABLES 


We consider the process schematically illustrated in 
Fig. m A pump beam at frequency uj p propagates 
through a quadratic nonlinear waveguide array. Nonlin¬ 
ear wave mixing combined with quantum fluctuations can 
convert a pump photon into two lower frequency photons 
called the signal and idler. The state of these photons is 
described by a biphoton wavefunction which evolves as 
they propagate through the array. The quantum corre¬ 
lations of the signal and idler photons leaving the array 
can be observed through coincidence measurements of a 
pair of single photon detectors. 

Theoretically, we employ the formalism of Refs. 0 , 111 , 
126], considering type -1 near-degenerate SPDC under con- 
tinuous wave pumping at frequency such that phase 
matching occurs when uj p ~ 2 c o s ,ii where uj s ~ uJi are 
the signal and idler photon frequencies. This can be im¬ 
plemented in experiments by placing an appropriately 
chosen spectral filter at the array output. 

Close to degeneracy, the waveguide coupling coeffi¬ 
cients for the signal and idler photons are approximately 
the same, ie. C s ,i = C. On the other hand, the higher 
frequency pump beam experiences much stronger con¬ 
finement and hence weaker coupling C p between neigh¬ 
bouring waveguides. A good approximation for recent ex¬ 
periments is C p ~ 0 , such that coupling of pump photons 
between waveguides can be neglected [19]. Consequently, 
under the undepleted pump approximation and choosing 
a frame rotating at the pump frequency, the pump beam 
profile remains constant along the waveguide array, sim¬ 
plifying the theoretical analysis considerably. 

Under these conditions, the evolution of the bipho¬ 
ton wavefunction can be described by a Hamiltonian 
H = + F( spdc ). jj(Q w ) accounts for the lin¬ 

ear diffraction (quantum walk) of biphotons through the 
waveguide array, and ( SPDC ) is a gain term accounting 


for their generation via SPDC. In normalized units with 
h = 1 , 


H (Q w) = £ [a^H mn a^ + H mn a$] , ( 1 ) 


ij(spDc) = h t + H. c ., (2) 

Tip 

where dn’^(dn’^) creates (destroys) a signal or idler 
photon at the nth waveguide, 7 is the nonlinear coeffi¬ 
cient, EnJ is the pump amplitude in waveguide n p , and 
H nrn are elements of the waveguide array’s tight bind¬ 
ing Hamiltonian. Diagonal elements H nn account for the 
propagation constant of the signal/idler photons in the 
nth waveguide; off-diagonal elements describe evanescent 
coupling between waveguides. 

We assume there is no decoherence or loss, such that in 
the absence of multiple photon pairs being generated si¬ 
multaneously, the biphoton state | \k) is pure and evolves 
according to the Schrodinger equation (l 2 l| . 


id z | *> 


tf( QW ) + tf(SPDC)j (| | } ( 3 ) 


where | 0) is the vacuum state. We will solve this equa¬ 
tion and reveal the effect of band structure topology by 
transforming to the eigenbasis of H^ w \ ie. the lattice’s 
Bloch wave basis. 

Consider a waveguide superlattice with a unit cell con¬ 
sisting of N waveguides. It is convenient to introduce the 
vector notation a n = (a n? 1 , ...,a n ,iv)> where now n num¬ 
bers the unit cell and a n ,m is the annihilation operator 
for the mth sublattice. Eq. © is recast as 


ff (QW| = E [<&>'».»* 


W 




i(0 


(4) 


where now each i7 mn is promoted to an N x N matrix, 
with off-diagonal elements accounting for coupling be¬ 
tween the different sublattices. Hmn = H m +\ n = H mn +i 
is periodic, such that transforming to reciprocal space 
a (k) = &nE kn puts into block diagonal form, 

H^)(k) = Y j H 0n e ikn . (5) 

n 

The eigenvectors of H^ w \k) are the superlattice’s 
Bloch functions u p (fc); eigenmodes of H^ w "> are Bloch 
waves, constructed by 

B p {k) = E ' u p( fc K fen > (6) 

n 

where • denotes the usual dot product, and p = 1, ...,7V 
is the band index. In this Bloch wave basis, FUQ W ) takes 
the simple diagonal form 

£(QW) = ^ T dk s (3 Ps (k s )Bl(k s )B Ps (k s ) 

Ps 

+ E r dkiPpi (h)Bp. (ki)B Pi (fcj), (7) 
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where /3 p (k) is the propagation constant of the Bloch 
wave in band p with crystal momentum k. To obtain 
i^(SPDC) [ n tpjg basis, we invert Eq. (j 6 ]), 

a n = ^2 f dku* p (k)B p (k)e~ lkn , ( 8 ) 

p 

and substitute into Eq. ©• Writing the pump ampli¬ 
tude in vector form in terms of its sublattice compo¬ 
nents, = (E nj i ,..., E n ^) : and applying a Fourier 

transform, E ^ p \k p ) = E ^e lkpn , we obtain 

fr(sPDC) = i j2 (dk s dkir p ,^(k a ,ki)Bl(k s )Bl(ki), 

Ps,Pi J 

(9) 

where 

N 

^p s ,pi ki) = 7 'y ] Ej (k s + ki)u Ps j (k s )u Pi j (&*), 
j =i 

(iO) 

is the coupling efficiency into the biphoton Bloch wave. 
The summation is over the N sublattices forming the 
superlattice. Eq. (J9]) is also diagonal in the Bloch wave 
basis, thus in a similar manner to Ref. [7] we can integrate 
Eq. ([3]) to obtain the output biphoton wavefunction (up 
to an overall normalization factor), 

!*> = £/ dk s dki T Ps . ( k s , ki ) Lsinc( A/3L / 2 ) 

p s ,Pi J 

x exp(-iA/3L/2)Bl s (k s )Bl i (k i ) | 0,0), (11) 

where A/3 = A/3® — /3 Ps (k s ) — /3 Pi (ki) is the phase mis¬ 
match into the biphoton Bloch wave, is the sin¬ 

gle waveguide phase mismatch, and L is the propagation 
length. 

Eq. (fill) tells us that two factors determine whether 
a biphoton Bloch wave is strongly populated: how close 
the mode is to phase matching (small |A/?|), and how 
strongly the pump beam profile is matched to the mode’s 
transverse profile (large |T|). Let us now discuss how the 
lattice topology can affect each of these. 

The phase matching condition A/3 = 0 depends only on 
the biphoton mode eigenvalues. Since the lattice topol¬ 
ogy is completely independent of the spatial dispersion 
/3(fc), in an infinite lattice the phase matching is insensi¬ 
tive to the topology: it cannot distinguish between two 
topologically distinct lattices 0 - On the other hand, in 
a finite lattice, nontrivial topology can result in “topolog¬ 
ically protected”, exponentially localized edge modes [ 2 l| . 
The phase matching condition is sensitive to these edge 
modes. 

The coupling efficiency |T| clearly depends on the 
Bloch function profiles u p (k) via Eq. (flQl) . Thus, we ex¬ 
pect nontrivial “winding” or topology of the Bloch func¬ 
tions to have some effect, even in an infinite lattice. 


While this Bloch wave decomposition is a convenient 
way to theoretically study SPDC in a superlattice, un¬ 
fortunately T and the Bloch functions u p (k) are not di¬ 
rectly observable in experiments. Instead, what is typ¬ 
ically measured is the magnitude of the biphoton wave- 
function Eq. ©ED, in either real or Fourier space, using 
coincidence measurements from a pair of single photon 
detectors. Therefore, instead of the Bloch wave basis in¬ 
dexed by band number p and crystal momentum k , we 
also need to consider the output in real and momentum 
space. For the latter, we use an extended Brillouin zone 
represenation, allowing k E [— Ntt, Ntt] to lie in the first 
N Brillouin zones, and using the Fourier amplitudes in 
the pth Brillouin zone as a proxy for the Bloch wave am¬ 
plitude in the pth band [27|, [28| . 

We would also like to quantify how “quantum” a given 
biphoton state is, and whether the lattice topology influ¬ 
ences the entanglement of the generated photons. One 
useful measurement of quantumness is the Schmidt num¬ 
ber [ 29 ], obtained via the singular value decomposition 

of |«), 

l = E vT l v j s) )® l T>’ ( 12 ) 
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where | Vj) are the Schmidt modes, A j > 0 are normalised 
such that JT Aj — 1 , and we define the Schmidt number 
as K — J2j 1/A|, which measures the number of entan¬ 
gled modes. 

III. TWO BAND MODELS 

The simplest case allowing for nontrivial topology is 
N = 2 band models, which often form a good approxi¬ 
mation to more complicated systems. The most general 
two band Bloch Hamiltonian Eq. (0 is 0 

= ( d, d +id u d ‘- d i dy ) - m ■ a, ( 13 ) 

where d = (a x , cr y ,cr z ) is a vector consisting of the three 
Pauli matrices, and we will see in the following that 
the vector d (k) = (d x (k),d y (k),d z (k)) provides a con¬ 
venient way to visualise both the Bloch functions and 
their topology. Diagonal elements of the Bloch Hamilto¬ 
nian H^ w \k) account for coupling between waveguides 
belonging to the same sublattice and their propagation 
constants, while off-diagonal elements account for cou¬ 
pling between different sublattices. 

Diagonalizing Eq. da, we obtain the Bloch wave 
eigenvalues and corresponding Bloch functions, 

P±(k) = ±|d(fc)|, (14) 

u ±(k) = sin ^1, ± cos ^ , (15) 

where we have introduced the spherical polar angles ip(k) 
and cj)(k) corresponding to the direction d (k) = d/|d| = 
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(a)Nontopological (b) Trivial (c) Nontrivial 



FIG. 2. (Color online) Bloch sphere representation of Bloch 
function topology. Bloch functions u p (k) form a closed loop 
on Bloch sphere surface as k traverses the ID Brillouin zone, 
eg. (a). (b,c) Topological phases phases are well-defined when 
u p (k) is constrained by symmetry to lie on a circle. Here chiral 
symmetry restricts u p (k) to the equatorial plane. The curve 
completely winds around the circle in the nontrivial phase (c), 
and it cannot be continuously deformed to (b). 

(cos (j> sin cp, sin </> sin cp, cos ip ). We assume the two bands 
are separated by a gap, ie. /3± ^ 0 , such that |d(fc)| ^ 0 
for all k and the angles </>, ip are always well-defined. 

The Bloch sphere provides a simple way to visualize the 
Bloch functions and their topology. The Bloch functions 
are mapped to points on the sphere’s surface specified by 
the pair of angles (</>, ip). Since the Bloch Hamiltonian is 
periodic, H^ w \k) = H^ w \k + 27t), <p(k) and p{k) are 
periodic as well. Therefore as k traverses the Brillouin 
zone, (</>, ip) maps out a closed curve on the sphere sur¬ 
face, see eg. Fig. ETa). This allows the topology of the 
Bloch functions to be mapped to the topology of curves 
on a sphere. 

The topology is trivial if the closed curve can be con¬ 
tinuously deformed to a point. In a general two-band 
model without any symmetries, such that d Xj y iZ are un¬ 
constrained and (</>, ip) can take any value, this is always 
possible and the topology is trivial. Hence nontrivial 
topology in a ID lattice requires some symmetry. An 
example is the “chiral” symmetry d z (k ) = 0 HH, which 
restricts (</>, ip) to the equatorial plane in Fig. [ 2 ](b,c) . As 
long as the two Bloch bands are separated by a gap, 
/3± 7 ^ 0 , there is no way to continuously deform the trivial 
curve (b) such that it has the nontrivial winding around 
the equator in (c). 

In a finite lattice, this nontrivial winding results in pro¬ 
tected edge modes with propagation constant /3 = 0 in 
the middle of the gap between the two bands. Their 
transverse profiles decay exponentially away from the 
edge of the lattice, at a rate determined by the size of the 
gap 0 . The modes are protected in the sense that they 
cannot be destroyed by any perturbation that respects 
the chiral symmetry d z = 0 as long as the two bands re¬ 
main separated by a gap 0 . The topological invariant 
associated with this protection is the Zak phase 0 - 

Turning to properties of biphoton modes, the signal 
and idler photons can excite different combinations of 
the two single photon bands with propagation constants 
given by Eq. nil. Thus there are four biphoton bands 
with energies ±(|d(fc s )| + |d(fc i )|) and ±(|d(fc s )| - |d(fci)|), 
corresponding to the signal and idler photons exciting the 


same and different single photon bands respectively. If 
signal and idler photons are near-degenerate, such that 
H s « Hi , then the \d(k s )\ - Jd(fc*)| and -|d(fc a )| + \d(ki)\ 
bands overlap, centred at /3 = 0. When the width of 
a single photon band is larger than the band gap, the 
biphoton spectrum is gapless. To show this, we note that 
A = 2min(|d(fc s )|) is the size of the single photon band 
gap and let w = max(|d(fc s )|) — min(|d(fc s )|) be the band 
width. Then the bottom of the top band, A, is below 
the top of the middle bands, (A/2 + w) — A /2 = w when 
w > A. 

There are two different types of biphoton edge modes, 
conventional and hybrid. In a conventional edge mode, 
both photons excite the same single photon edge mode; 
thus /3 = 0. This mode is inevitably degenerate with 
modes belonging to the middle pair of biphoton Bloch 
bands, which are centred at /3 = 0. 

In a hybrid mode, one photon excites the edge mode, 
while the other excites a Bloch band mode; thus the 
hybrid modes form bands with /3 = ±|d(fc s )| and /3 = 
±|d(fci)|. When the biphoton Bloch wave spectrum is 
gapped, hybrid modes with w < |d(£^)| < A reside in 
the gap. 

We obtain the coupling efficiency by substituting the 
mode profiles Eq. (fl5l) into Eq. m , 

r p . tPi (k e , hi) = E[ p \k s + ki) sin(^/ 2 ) sm^ i /2)e~ i ^ + ^ 

+ p s PiEff\k s + ki) cos(<£> 5 /2) cos(<pi/2), 

(16) 

with p Sri = ±1 and p s ,iAs,i = <p(k Si i), </>(k 8ii ). Recall 
E[ p (fc) is the Fourier transform of the pump amplitude 
on the two different sublattices. p s pi = 1 (— 1 ) if signal 
and idler photons come from the same (different) Bloch 
bands. This sign can always be absorbed into the relative 
phase of e[ p I , which means that, as far as the coupling 
efficiency is concerned, all the biphoton bands look the 
same, so the only difference will be in their dispersion 
(phase matching conditions). 

Defining the vector 

n = (sm((p s /2) sin((^/2)e _ ^ s+< ^, 

PsPiE!f\k s + ki) cos(ip s /2) cos(<pj/2)) , (17) 

the angles {ip s ,i,(j)s,i) define a direction on the Bloch 
sphere, and Eq. m can be recast as T = E ( p \k s + 
ki) -n. Hence, the coupling efficiency is maximized when 
E W(k 3 + k i) is parallel to n, and zero if it is perpendic¬ 
ular. 

In a trivial phase p s y and (j) s ^ do not exhibit any wind¬ 
ing. If the gap is sufficiently large, they typically stay 
close to mean values independent of k 8j i, eg. Fig. [21(b). 
Then it is possible to shape the pump profile E ( p \k s ,ki) 
such that all the modes in a band are strongly excited, 
resulting in behaviour similar to the homogeneous lat¬ 
tice case 0. On the other hand, in the nontrivial phase 
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the winding of (p 3j i or (f) s ^ means that it is impossible to 
simultaneously excite all modes efficiently: shaping the 
pump profile such that is parallel to n for some k s ^, 
there are inevitably other values of k s ^ for which they are 
perpendicular. Thus, there are selection rules preventing 
the excitation of some modes. This is the main conse¬ 
quence of nontrivial winding or topology in an infinite 
lattice. 

So in summary, in ID two band lattices the nontrivial 
topology has two main effects: 

• when the single photon band gap is sufficiently 
large, there are hybrid biphoton edge modes with 
frequencies lying in the band gaps of the Bloch wave 
spectrum 

• the coupling efficiency for the excitation of Bloch 
waves is modulated, giving additional selection 
rules for the generation of biphotons 

In the next Section we apply these ideas to a concrete 
example. 


IV. SU-SCHRIEFFER-HEEGER MODEL 


The Su-Schrieffer-Heeger (SSH) model [HI, 32] presents 
a simple example of a ID topological phase. Physically, 
it describes a ID waveguide array where the waveguide 
separation is modulated such that the nearest neighbour 
coupling strength alternates between C+SC and C — SC , 
see Fig.[3](a). This is described by the single photon tight 
binding Hamiltonian 

ff (QW) = £ ([C+{-l) n 5C\ai +1 a n 

n 

+[C-{-l) n 5C]ai_ l a n ), (18) 

Introducing two sublattices and applying a Fourier trans¬ 
form, we obtain the single photon Bloch Hamiltonian 

H^ w \k) = [C-\-SC-\-(C—SC) cos k\cr x +(C—SC) smkd y , 

(19) 

which corresponds to Eq. m with d z = 0. The spec¬ 
trum, 


p±(k) = ±V2 ye 2 + W 2 + Jc + 8C){C - SC) cos jfe, 

( 20 ) 

is plotted in Fig.[3](b) in the extended Brillouin zone rep¬ 
resentation. When SC = 0 the model reduces to a ho¬ 
mogeneous lattice with the usual single dispersion band. 
Nonzero SC doubles the lattice period, forming two sub¬ 
lattices and splitting this band in two, each with a width 
of 2 (|C| — \SC\), and separated by a gap of size 4\SC\. 

Single photon Bloch functions are obtained as the 
eigenvectors of Eq. 


U ±(k) 


1 

71 

i 

71 


(l,±l[(C-dC)e~ ik + (C + 5C)}) , 
(e _ ^ (fe) ,±l) , ( 21 ) 



FIG. 3. (Color online) The ID SSH model, (a) Lattice geom¬ 
etry with modulated waveguide separation (coupling), unit 
cell marked by dashed line, (b) Spectrum, shown in extended 
Brillouin zone scheme, for homogeneous (SC = 0, dashed) 
and dimerized (\SC\ = C /2, solid) lattices, (c) Phase wind¬ 
ing of eigenmodes in trivial (SC = C/2, blue) and nontrivial 
(SC = —C/2, red) phases. 


so recalling Eq. ([T5]h ip(k) = 7 r/ 2 , ie. Bloch functions live 
on the Bloch sphere’s equatorial plane. The winding of 
4>(k), plotted in Fig. [3](c) , determines the lattice topology. 
In the trivial phase SC > 0 , intracell coupling is stronger 
and <p(k) shows no winding, with cj)(k) ~ 0 for all k. In 
the nontrivial phase SC < 0 , intercell coupling is stronger 
and cj)(k)—(j)(k-\- 27t) = 27r, winding once around the Bloch 
sphere’s equatorial plane. Note that in the limit of an in¬ 
finite lattice, the choice of unit cell boundary is arbitrary, 
so this “topology” becomes ill-defined. However, in any 
finite lattice, the nontrivial phase is distinguished by a 
pair of “topologically protected” edge modes at /3 = 0 
(the middle of the band gap) [H]- These modes appear 
when a strong bond is broken to form the edge. Experi¬ 
mentally, one can compare the two phases using a single 
lattice with the two edges terminated differently. 

Biphoton modes /3 ns?n . (& s , &i) are constructed as com¬ 
binations of pairs of single photon modes. The biphoton 
spectrum is shown in Fig. 0Ja), once again using the ex¬ 
tended Brillouin zone representation. In contrast to the 
single photon case, here the spectrum remains gapless for 
\SC/C\ < 1/3 (as long as the width of a single photon 
band exceeds the band gap). 

The nontrivial phase hosts pairs of biphoton bands 
at each edge corresponding to the “hybrid” edge modes, 
shown in Fig. 0|(b). These bands bifurcate from the the 
topological phase transition at SC = 0 . However, for 
—1/3 < SC/C < 0 they still overlap with the bulk Bloch 
bands. Notice how their edges remain pinned at /3 = ± 2 . 
This is because the “topological protection” ensures the 
single photon end mode remains fixed at /3 = 0. In con¬ 
trast, the trivial phase SC > 0 does not have any edge 
modes 


A. Pumping - infinite lattice 

Here we consider in-phase pumping of two adjacent 
waveguides far from the lattice edges (the infinite lat¬ 
tice limit). We will start by considering examples of nu¬ 
merical solutions of Eq. m before discussing how the 
results generalise. In the following examples, we use 
\SC/C\ = 1 / 2 , and assume a normalised propagation dis¬ 
tance L = 10/C. 
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(b) Density of states 

4. 



0 SCIC 1 


max 
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FIG. 4. (Color online) (a) Biphoton Bloch wave spectrum, 
shown in extended Brillouin zone scheme, SC/C — 1/2. (b) 
Density of states as a function of coupling modulation 5C/C. 
Bloch bands (solid red region) are gapped for \8C/C\ > 1/3. 
Edge modes only exist in the nontrivial phase, SC/C < 0. 


Trivial 

pump 


JldJULTUTL 



Nontrivial 

pump 


numiin 



-271 £ 271 

s 


FIG. 5. (Color online) Real and Fourier space biphoton cor¬ 
relations, trivial (left) and nontrivial (right) phases. \SC\ = 
C/2. LC — 10. A/3 = 3 C (band 1 excitation). 


Fig. [5] shows the output biphoton correlations in the 
trivial and nontrivial phases when the pump frequency 
is resonant with the first biphoton band, A/3^ = 3(7. 
The output in the trivial phase displays bunching and 
antibunching, resembling the output of a homogeneous 
lattice with a single waveguide pump [7j. In contrast, 
pronounced antibunching occurs in the nontrivial phase; 
photon bunching is strongly suppressed. 

This result can be understood quite intuitively by con¬ 
sidering the strong modulation limit \SC\ « ( 7 , in which 
the lattice consists of strongly coupled “dimer” pairs of 
waveguides, with weak coupling between neighbouring 
dimers. At each dimer, the single waveguide modes 
hybridize to form in- and out-of-phase modes. When 


Trivial Nontrivial 

pump pump 


nrjuLTLn rm_rm_njL 



FIG. 6. (Color online) Real and Fourier space correlations, 
trivial (left) and nontrivial (right) phases. \SC\ = (7/2. LC = 
10. A/3 = 0 (bands 2 &3 excitation). 


SC > 0, the pump excites the in-phase mode of a single 
dimer; thus the output resembles that of a homogeneous 
lattice when a single waveguide is excited. In the non¬ 
trivial phase, the pump excites two dimers; interference 
between them suppresses photon bunching. 

The response changes when the pump is tuned to the 
2nd and 3rd (overlapping) bands. The correlations in 
Fig. [6] no longer show any significant qualitative differ¬ 
ence between the trivial and nontrivial phases; both dis¬ 
play strong antibunching. The main quantitative differ¬ 
ence is that the total output intensity is an order of mag¬ 
nitude smaller in the trivial phase. 

Let us now relate these observations back to the gen¬ 
eral theory. The pump detuning A/3 imposes a phase¬ 
matching condition on the spatial modes, such that only 
Bloch modes in resonance can be strongly excited. How¬ 
ever, the Bloch mode spectra for SC and —SC are iden¬ 
tical, so the coupling efficiency T is solely responsible for 
the differences in the biphoton correlations. Evaluating 
Eq. ([T6]b we obtain 

r = I (F P) (fc s + fcJe-W^+W*)] + PsPi E {p \k a + kij) , 

( 22 ) 

where p s ^ = d=, and \ (k s + ki) is the Fourier transform 
of the pump amplitudes on the two sublattices. Since 
only a single unit cell is pumped, the Fourier transform 
is a constant, ie. E^i/kg+ki) = e[ p I , independent of k s ^. 
We see that the topology, via the phase (j){k s5 ^), affects 
the interference between the two sublattices, which in 
turn controls T. When a single sublattice is pumped (eg. 
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(a) Nontrivial (b) Trivial (c) Nontrivial Trivial 



FIG. 7. (Color online) SPDC coupling efficiency |r(fc s , ki)\ for 
two waveguide pumping, normalised by pump intensity, equal 
pump intensities, 7 r relative phase between the two waveg¬ 
uides, in nontrivial (a) and trivial (b) phases. \SC/C\ — 1/2. 
The nontrivial phase displays a maximum, and minimum of 
zero along curves in 2D Brillouin zone. In trivial phase, these 
lines are not present, |F| is controlled mainly by the rela¬ 
tive pump phases, (c) Contrast in the coupling efficiency, 
max(|r(fc s , ki)\) - min(|r(fc s , fc)l), as a function of relative 
pump phase 0 and coupling modulation SC 


E^ = 0), no interference occurs, so \T\ is independent 
of </> and SC and the topology is irrelevant. So both 
sublattices must be pumped to observe any sensitivity to 
the topology. 

In the nontrivial phase, as either k s or ki traverse the 
Brillouin zone, 4>(k S: i) takes all possible values. Hence 
there is always a curve through the (k s ,ki) Brillouin 
zone where |T| attains its maximum of ( \e[ p ^ | + {E^ |)/2 
(constructive interference), and one where |F| drops to 
its minimum of (\e[ p ^\ — (destructive interfer¬ 

ence). When \e[ p ^\ = \E^\ this minimum is zero, so we 
have perfect destructive interference and coupling into 
the corresponding Bloch wave is forbidden. Changing 
the relative phases of e[ p shifts the two curves, but it 
does not remove them. In contrast, in the trivial phase 
<f>(k Si i) does not display any winding, 4>(k S: i) « 0, and the 
efficiency is only controlled through the relative phase 6 

of e[% 

We demonstrate these two different cases by plotting 
\T(k s , ki)\ in Fig. [7|[a, b), assuming out-of-phase pumping 
of the two waveguides and phase matching with band 
1. In this case, |F| vanishes for both phases when ki = 
— k s , ie. antibunching is suppressed. In the trivial phase 
however, |F| remains small for all (, k s ,ki ), ie. no modes 
are efficiently excited, while |T| attains its maximum of 1 
in the nontrivial phase and efficiently excites some modes. 
This explains the results in Fig. [6] (similar correlations, 
but different intensities). 

More generally, we show in Fig. [71(c) the effect of 
SC and 0 by plotting the contrast, max(|r(& s , ki)\) - 
min(|r(fc s , ki)\), for the case \e[ p ^\ = |E^|- We verify 
the reasoning in the previous section that in the nontriv¬ 
ial phase, the contrast is always maximum, while in the 
trivial phase it decreases to zero. 

The total contrast controls the efficiency of the SPDC. 
Fig. [51(a) shows the total down-converted intensity as a 
function of SC/C and 6. In the nontrivial phase, there 


(&) Total intensity (b) Schmidt number 



0 5C/C 1 -1 0 5 c/C 1 


FIG. 8. (Color online) (a) Total intensity of down-converted 
photons and (b) Schmidt number as a function of SC/C and 
relative phase 0 of the two pumped waveguides. A/3^ = 
2(|C| + 15(71) (middle of band 1 excitation). 


are always some modes that are spatially matched with 
the pump and thus strongly excited. Hence the total out¬ 
put intensity is relatively insensitive to the relative pump 
phase 0. In the trivial phase, it is crucial that 6 matches 
the selected band’s Bloch wave profile, otherwise no spa¬ 
tial modes are strongly excited and the down-converted 
intensity vanishes. 

Fig. [51(b) shows the influence of 6 and SC/C on the 
Schmidt number of the biphoton state, assuming the 
pump frequency is tuned to the centre of band 1. The lat¬ 
tice topology has a less significant effect here; in fact the 
Schmidt number is largest when SC « 0 (a homogeneous 
lattice). 

So far we have focused on a pump that is confined to a 
single unit cell of the lattice (two waveguides). Let us now 
consider briefly the effect of a broader pump beam. As 
the pump profile is made wider, it becomes more localized 
within the lattice’s Brillouin zone, such that | E^f\ | = 

| e[ p I I (k s + ki) acquires a ^-dependence. 

Applying Eq. ([22lb this localization of the pump 
around some point in the Brillouin zone leads to an addi¬ 
tional, k Sj i -dependent modulation of |F|. Notice however 
that this modulation is distinct from that arising from 
the Bloch waves themselves: it depends on the sum of 
the down-converted photon wavenumbers, k s + ki, in¬ 
stead of k Sr i individually. Furthermore, this additional 
modulation induced by a broad beam is independent of 
the lattice properties, including its topology. 

As an example, we consider two examples of bipho¬ 
ton correlations when two unit cells (four waveguides) 
are pumped. Pumping the waveguides in phase in Fig. [9] 
favours antibunching [7J. The selection rule imposed by 
the lattice topology (which favours antibunching only 
in the nontrivial phase) becomes redundant, and both 
phases display similar correlations and antibunching. 
Conversely, pumping the two unit cells with a 7 r phase dif¬ 
ference in Fig. [TOl promotes photon bunching [7|], clearly 
visible in the trivial phase. In the nontrivial phase, this 
additional selection rule competes with topology-imposed 
suppression of bunching to produce a complex pattern of 










8 


Trivial Nontrivial 

pump pump 
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FIG. 9. (Color online) Real and Fourier space correlations, 
trivial (left) and nontrivial (right) phases. \SC\ = C/2. 
LC = 10. A/3 ^ — 3 C (band 1 excitation). Four waveguides 
pumped. 


Trivial Nontrivial 

pump pump 
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FIG. 10. (Color online) Real and Fourier space correlations, 
trivial (left) and nontrivial (right) phases. \SC\ = C/2. LC — 
10. A/3^ = 0 (bands 2&3 excitation). Four waveguides 

pumped, 7r phase difference between unit cells. 


correlations. 


B. Pumping - edge of lattice 

We next consider a finite lattice and the role of the 
topologically protected edge modes. Fig. fTTta-d) shows 
the real space output biphoton correlations when the 
waveguide at the end of the lattice is pumped, for differ¬ 
ent modulation strengths SC and pump detunings 

First, when the pump is tuned to the conventional edge 
mode, A/3^ = 0, we observe strong localization of the 
output in Fig. [TTT a). even though the pump is also res¬ 
onant with Bloch waves. This is because of the strong 
overlap of the pump profile with the edge mode. 

If the pump is tuned to the centre of the hybrid edge 
mode band, A f3^ = yj2(C 2 + SC 2 ), there are four dis¬ 
tinct regimes depending on SC: 

• SC/C > 1/3 (trivial and gapped): The pump is 
tuned to a band gap, so no modes are resonantly 
excited. 

• 0 < SC/C < 1/3, (trivial and gapless). The pump 
resonantly excites bulk biphoton modes, which 
propagate away from the edge. 

• —1/3 < SC/C < 0, (nontrivial and gapless). The 
pump resonantly excites bulk modes and an edge 
mode. 


• SC < —1/3, (nontrivial and gapped). Only a 
biphoton edge mode is resonantly excited. One 
photon is trapped at the edge, while the other prop¬ 
agates into the bulk. 

Fig. fTTl b.c) demonstrates the last two regimes. 

For comparison Fig. fTTl d) shows the output correla¬ 
tions when SC = 0 (homogeneous lattice) and there are 
no edge modes. Photon bunching occurs, with signal and 
idler both propagating into the bulk. 

We consider more generally in Fig. HHe,f) how the 
biphoton intensity and Schmidt number depend on the 
pump detuning and coupling modulation. Due to the 
strong overlap with the pump beam, the output inten¬ 
sity is maximum when the pump is resonant with the 
conventional edge mode. However, since only a single 
mode is strongly excited, the Schmidt number K « 1 re¬ 
veals there is no entanglement. Similar to the bulk case, 
K is maximized for relatively small SC, when the pump 
is tuned to the centre of the Bloch bands. 

In summary, we have shown how the SSH model can 
exhibit nontrivial topology in biphoton correlations: in 
the bulk (additional selection rules), and at the edge (“hy¬ 
brid” biphoton edge modes). 


V. MODULATED LATTICE DEPTH 

For comparison, we briefly consider here an experimen¬ 
tally accessible “nontopological” model: a binary lattice 
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FIG. 11. (Color online) (a-d) Real space biphoton correlations 
when edge of lattice is pumped. Parameters indicated in (e). 
(a) Pump resonant with conventional edge mode, (b) Pump 
resonant with hybrid edge mode, (c) Pump resonant with 
hybrid edge mode and Bloch waves, (d) Homogeneous lattice 
with no edge modes, (e) Total biphoton intensity and (f) 
Schmidt number of the biphoton state as a function of pump 
detuning and coupling modulation. 


where the waveguide depths are modulated with strength 
m, while the waveguide spacing and coupling strength C' 
are constant, see Fig. [T2l 

Following the same procedure as for the SSH model, 
we obtain the Bloch Hamiltonian 

iI (QW) (fc) = C"(1 + cos k)a x + C' sin kd y + ma z , (23) 
with spectrum 

0± = ± v /2C ,,2 (l + cosfc) +m 2 , (24) 

and Bloch functions 

u ±(*0 = (m±/3,C'(l + e ik )), (25) 

N 2 {k) = m(m ± (3) + 2C /2 (1 + cos k). (26) 

We set m = 2 SC and C' = VC 2 — SC 2 , so that the bulk 
spectrum is identical to the SSH model’s Eq. (l20lb allow¬ 
ing for a fair comparison. The only remaining difference 


—777 +m 

ju uJ 


uni l 



FIG. 12. (Color online) Biphoton correlations in a binary 
lattice, tuned to have same Bloch wave spectrum as the SSH 
model with 5C/C = 1 /2. (a,b) A = 3 C, (c,d) A 0^ = 0. 
Compare against Fig. [5] and Fig. [6] respectively and observe 
the resemblance to the trivial phase. 


between the two models is the singular winding of the 
Bloch functions in the SSH model, which is absent here. 

At k = 0, the Bloch functions excite both sublattices, 
while at k = n (Brillouin zone edge), they reside on a 
single sublattice only, with energies ±m. Fig. [21(a) shows 
the Bloch functions using the Bloch sphere representa¬ 
tion. Since the curve can be continuously shrunk to a 
point, there is no nontrivial winding as the Brillouin zone 
is traversed. This is because taking the limit m oc 
continuously deforms the lattice to an effectively homo¬ 
geneous chain without closing the band gap: one can see 
this by Taylor expanding the dispersion relation Eq. (|24| ) 
for small C jm and recovering a /3 ~ cos k dispersion re¬ 
lation, with the Bloch functions independent of k. In a 
finite lattice, there are no edge modes. 

We consider similar to Figs. 15161 pumping two adjacent 
waveguides far from the lattice edge. Similar to the triv¬ 
ial phase in the SSH models, the output correlations in 
Fig. [T2l resemble those from pumping a single waveguide 
in a homogeneous lattice: pumping the 1st band reveals 
both bunching and antibunching, while antibunching is 
favoured in bands 2 and 3. 

This example further highlights how lattices can pro¬ 
duce very different biphoton correlations even when their 
spectra (eigenvalues) are identical and they are pumped 
in the same way, because of the additional selection rules 
imposed by the Bloch functions and their trivial or non¬ 
trivial topology. 
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VI. CONCLUSIONS AND OUTLOOK 

In summary, we have explored the effect of lattice 
topology on spontaneous parametric down-conversion in 
one dimensional quadratic nonlinear waveguide arrays. 
We have shown how nontrivial winding in the Bloch wave 
spectrum leads to selection rules for the generation of 
entangled photon pairs. Finite lattices can host topo¬ 
logically protected edge modes, which interestingly en¬ 
able the generation of entanglement between localized 
and propagating spatial modes. As a specific example 
we considered in detail an analogue of the Su-Schrieffer- 
Heeger model, which can be experimentally realised in 
lithium niobate nonlinear waveguide arrays using exist¬ 
ing fabrication techniques. 

The study of extensions to two dimensional topolog¬ 
ical phases remains an open problem. Can the bipho¬ 
ton spectrum /3(k s ,ki ) and its eigenmodes host genuinely 
two dimensional effects, such as nonzero Chern number? 


Presumably, the “edge modes” in such a system would 
involve one photon bound at the edge, with the other 
propagating into the bulk. Another possible avenue to 
explore is SPDC in two dimensional waveguide arrays 
with nonzero Chern number. A two dimensional array 
results in a four dimensional biphoton spectrum, which 
raises the intriguing possibility of emulating highly exotic 
topological phases, such as the four dimensional quan¬ 
tum Hall effect [33|. However, experimental realisations 
would be quite challenging, since so far SPDC has been 
limited to one dimensional nonlinear waveguide arrays. 
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